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A MonteCarlo simulations study has been performed in or- 
der to study the effect of Pb as surfactant on the initial growth 
stage of Co/Cu(lll). The main characteristics of Co grow- 
ing over Cu(lll) face, i.e. the decorated double layer steps, 
the multiple layer islands and the pools of vacancies, disap- 
pear with the pre-evaporation of a Pb monolayer. Through 
MC simulations, a full picture of these complex processes is 
obtained. Co quickly diffuses through the Pb monolayer ex- 
changing place with Cu atoms at the substrate. The exchange 
process diffusion inhibits the formation of pure Co islands, 
reducing the surface stress and then the formation of multi- 
layer islands and the pools of vacancies. On the other hand, 
the random exchange also suppress the nucleation preferential 
sites generated by Co atoms at Cu steps, responsible of the 
step decoration. 



PACS numbers: 68.35.Fx, 68.55.-a, 8F10.-h 



I. INTRODUCTION 

The startling, and sometimes surprising physical prop- 
erties observed in the new series of artificial materials, has 
catapulted an intensive work around the world. The ex- 
perimental conditions needed to obtain atomically clean 
and nearly perfect flat substrates are achieved in most 
surface labs. Currently, almost any atom can be de- 
posited over almost any surface, allowing the creation 
of phases that are not found in nature. But, the way the 
atoms diffuse over the terrace (intra layer diffusion) and 
the ability to overcome terrace steps (inter layer trans- 
port) before binding other walking atoms, will determine 
the morphology of the growing film. Being intermixing 
and surface roughness development the most common ar- 
tifacts in the preparation of these materials, the obtaition 
of flat films through layer by layer growth (LbL) is mostly 
the exception than the rule. External parameters, like 
substrate temperature and evaporation rate among the 
most usual, and internal ones, like interaction potentials, 
Ehrlich-Schwoebel barrier, diffusion type, surface orien- 
tation and so on, are crucial in determining the type of 
film growth |l|-Ql . On the other hand, the evaporation or 



adsorption of some materials, known as surfactants, pre- 
vious to the metal deposition has proved to promote LbL 
growth A perfect example of this behavior is the 

growth of Cu over Cu(lll). In fact, while in absence of 
surfactant Cu grows in a 3D fashion over Cu(lll), after 
the pre evaporation of a monolayer of Pb, LbL growth 
in a wide range of temperatures and evaporation rates 
can be obtained, [gj In a recent work, we showed that 
the monolayer of Pb changes the diffusion mechanism of 
Cu over the surface from hopping to atomic concerted 
exchange, turning the diffusion over Cu(lll) surface like 
over the Cu(100) one, where LbL growth is usually ob- 
tained @. 

The growth of Co over Cu(lll), on the other hand, is a 
good example of hetero-epitaxial non LbL growth. Mul- 
tiple height island formation, step decoration, vacancy 
(pools) islands and Cu-Co intermixing are common fea- 
tures in this system [p|JTc|]. In a recent work, we showed 
that this behavior depends on the adsorbate - substrate 
interaction, and it should be awaited any time materials 
of these types are put into contact . As in the case of 
Cu/Cu(lll), the pre-evaporation of a monolayer of Pb 
acts as a surfactant promoting 2D growth, and allow- 
ing the fabrication of good quality magnetic multilayer. 
In this work we present results, based on MonteCarlo 
simulation about the effect of Pb as surfactant for het- 
eroepitaxial growth, in particular we analyze the case 
Co/Cu(lll) 



II. COMPUTATIONAL DETAILS 

The interaction potentials employed in our MC simu- 
lations have been obtained using the second-moment ap- 
proximation of the tight-binding scheme (TB-SMA) p2| . 
They include a short-range, repulsive pair potential plus 
a long-range, many-body contribution based on a tight- 
binding description of the electronic structure |]l3|l . They 
have been successfully used previously to study the prop- 
erties of transition metal surfaces jljjl^], noble metals 
|l6| and microclusters ||l2] , [l7| , |l8| . This potential will be 
used for the six involved interactions, i.e. Co-Co, Pb-Pb, 
Cu-Cu, Co-Cu, Co-Pb and Cu-Pb with their correspond- 



ing physical parameters. JTT] ] The potential at the i — th 
atom is given by 

Vi(nj) = U{A^2exp[-p(nj - r )] 
i 

-[^eM-Mnj-ro)]] 1 / 2 } (1) 



where is the distance between the two atoms at n 
and rj, respectively, tq is the nn distance of the bulk 
crystal, U, A, p and q are parameters. In this paper, 
the sums are performed over all atoms within a sphere of 
radius equal to 3.2 times the Cu nn distance. We have 
checked that including interaction between further atoms 
does not alter the conclusion of this paper. The many- 
body interactions are included in the potential through 
the square root of the second sum. This term takes 
into account the essential band character of the metallic 
bonds The value of A is determined mini- 

mizing the cohesive energy of the bulk crystal with nn 
distance ro, while the values of U, p and q are deter- 
mined in such a way that the bulk cohesive energy and 
the bulk modulus calculated using (1) and the experi- 
mental value of ro are in good agreement with the ex- 
perimental bulk cohesive energy and bulk modulus. The 
values p — 9/vq and q — 3/ro are found to be appro- 
priate for the surface jlij . The interactions between the 
different types of atoms are determined by sets of six pa- 
rameters {Ui, pi, qi}; i = 1 stands for Cu-Pb pairs, i = 2 
for Pb-Pb, and so on. r$ i are the bulk nearest-neighbor 
(nn) distances. The validity of this approximation has 
been tested in previous studies with this combination of 
elements. The ordered structures and melting behav- 
ior of Pb overlayers on Cu(100) have been successfully 
reproduced In recent studies, the bulk properties 

of Pb, including the phonon spectrum, have been calcu- 
lated giving excellent agreement with the experimental 
results |2^ , |23|| . The Cu(lll) sample used in our calcu- 
lations consisted of 6 layers with 128 Cu atoms in each; 
to study the surfactant effect, a compact layer of 72 Pb 
atoms was added to the upper surface of the slab. The 
two lowest Cu layers were frozen to simulate the bulk, 
while all the remaining atoms in the sample were allowed 
to move. Most of the calculations have been performed 
at a sample temperature of 540 K in order to speed up 
the diffusion processes and improve the statistics. 



III. RESULTS AND DISCUSSION 

As different as the growth processes of Cu and Co over 
Cu(lll) is, the effect of Pb as surfactant recognizes the 
same root. Cu grows over Cu(lll) in 3D from the begin- 
ning. This is so due to the low intralayer diffusion barrier 
compared to the interlayer one (Ehrlich-Schwoebel bar- 
rier). The presence of a monolayer of Pb changes the 



diffusion mechanism from hopping to exchange, making 
both barriers of the same order. In absence of any other 
effect, this is enough to promote 2D growth M. The het- 
eroepitaxial growth of Co over Cu(lll) is by far more 
complex. The growth of Co over nude Cu presents the 
formation of double atomic height decorated steps, mul- 
tilayer islands and pools of vacancies @. Through ion 
scattering spectroscopy |24 and chemical titration pH] 
it has been shown that these islands are actually formed 
by a mixture of Co and Cu atoms. In a recent paper 
|l0f we showed that the diffusion mechanisms, the larger 
surface energy of Co as compared to Cu, and the Co-Co 
and Co-Cu interaction potentials are the responsible of 
these features. Over the Cu terrace Co diffuses through 
hopping. In this way, the adatom nucleation gives place 
to the formation of pure Co island. Once the Co islands 
reach a critical size, they explode forming multilayer al- 
loyed islands, and leaving pools of vacancies. On the 
other hand, the interlayer mass transport is produced, as 
in the case of Cu-Cu, through the exchange mechanism. 
But, in this case the process gives place to alloyed steps 
that act as nucleation sites for the diffusing Co adatoms, 
forming the double atomic height decorated steps. 

The question to answer is then if the effect of Pb 
on Co/Cu(lll) growth is similar to the homoepitaxial 
case. And if the change in the diffusion mechanism is 
enough to suppress all the different features appearing 
in Co/Cu(lll), promoting 2D growth. This is just the 
case; in Fig.l we show the evolution of a Cu (left panel) 
and a Co adatom (right one) over a Cu(lll) surface pre 
- covered with a Pb monolayer, obtained through MC 
simulations. In the upper panel the evolution of the z 
coordinate (normal to the surface) is depicted vs. MC 
steps. As we know, a drop in the Z coordinate is a fin- 
gerprint of a concerted exchange process. Thus, the first 
drop observed in both cases represents the exchange of 
Cu (Co) atom with the surfactant. This exchange is very 
fast and it is the responsible of the floating of the Pb 
monolayer. The second drop in z coordinate (not always 
occurring within the elapsed time analyzed) corresponds 
to the CoCu (CuCu) exchange. As we showed for Cu 
S, the Co diffusion over Pb/Cu(lll) system occurs be- 
low the Pb layer and through the exchange mechanism. 
We should note that although Co and Cu adatom dif- 
fusion below Pb are quite similar, they are not identical 
as well. In the lower panel of Fig.l we show a couple 
of snapshots of the MC diffusion simulations. We can 
see that substrate atoms involved in the concerted ex- 
change mechanism are no the same in both cases. While 
CuCu exchange involves next nearest neighbors, the pro- 
cess in the CoCu case is between nearest neighbors. This 
difference is maintained for the following steps, always 
in the vicinity of the implanted Co atoms. Thus, al- 
though the adatom diffusion below the surfactant layer 
is still through concerted exchange, the presence of Co 
as impurity modifies the nature of the exchanging mech- 



anism. This difference with Cu/Cu exchange may have 
important physical consequences. In fact, we have al- 
ready shown [jlO| that the interaction CoCu produces a 
dip binding well, locating the evolving atom in the near 
defect zone. Thus, the hopping mechanism will be more 
restricted than in the case of Cu/Cu. In fact, we are re- 
sembling the situation responsible of the step decoration 
p0| , but in the middle of the terrace. Additionally, the 
probability of a Cu Co reverse exchange, i.e. extracting 
the Co atom again, will be larger than in a "long jump 
exchange "like in CuCu case, increasing in this way the 
probability of reacher Co island formation. 

The change in the diffusion mechanism, from hopping 
to exchange is, as in the case of Cu over Cu(lll), enough 
to explain all the features of the surfactant effect over the 
heteroepitaxial growth. In fact, through the exchange 
mechanism, very quickly the diffusing atoms turn to be 
Cu ones, and no pure Co island could be formed. The al- 
loyed islands lowers the surface energy affecting the crit- 
ical size instabilities we found for Co over nude Cu ]To| ]. 
The non-adiabatic production of these alloyed islands, 
coming from the unstable pure Co islands, was the re- 
sponsible of the appearance of the pools of vacancies and 
the multiple height island, which are hardly formed in 
the new scenario. The supression of the decorated steps 
is based on the same mechanism. Since the diffusing 
atoms are mostly Cu ones, the CoCu alloyed steps, and 
thus the preferential nucleation sites located there, are 
not formed. The interlayer diffusion is still through con- 
certed exchange, but the morphology of the steps does 
not change, since the exchanging process is mostly CuCu. 

We have established that the 1 surface energy is the 
reason for the explosion of the Co islands., i.e. pure Co 
islands over Cu(lll) surface are naturally unstable. The 
exchange diffusion mechanism will tend, as we discussed 
in the previous paragraph, to the formation of mixed Co- 
Cu islands, lowering the surface energy. But, we have not 
shown still to what extent this effect is sufficient to sup- 
press the double height island formation. In fact, it could 
occur that the double height islands are still formed un- 
der different experimental conditions. For instance they 
may form with larger critical size, or they could depend 
on either coverage or evaporation rate. Actually, we have 
found that the CoCu islands are stable below the surfac- 
tant layer in all cases, even for the pure Co islands, due 
to the stabilizing effect of the Pb layer. In Fig. 2 we show 
this effect through a couple of snapshots of the evolution 
of a Co island contains 28 atoms with and without the 
Pb covering layer (Pb atoms are not shown for clearness). 
The first (two upper) panel show the initial situation, the 
same for the with and without surfactant case. The size 
of the studied island is quite larger the critical one jIo| . 
In the second panel, corresponding to Co over Cu(lll) ( 
middle) the complete process of multilayer mixed island 
and the appearance of vacancy pools is apparent. In the 
lateral sight (the right one) the mixed composition of the 



island is evident. In the upper vision (left one), the mo- 
bility of the island is also apparent. In the third lower 
panel, the surfactant assisted case, the Pb stabilized is- 
land is shown. Thus, the Pb modification of the diffusion 
mechanism not only prevent the formation of pure Co is- 
lands, promoting 2D growth, but also, in the remote, but 
statistically possible case that either Co pure or Co rich 
island are formed, the Pb layer acts stabilizing them. 

IV. CONCLUSIONS 

Through MC simulation, we have studied the growth 
of Co over Cu(lll) in the presence of a surfactant (Pb) 
layer. We found that the surfactant changes the diffu- 
sion mechanism of Co over the surface from hopping to 
concerted exchange. An important difference with the 
previously reported surfactant effect on Cu/Cu growth, 
is that the exchange process involves nearest neighbor 
atoms. This effect reduces even more the hopping mech- 
anism, and increase the probability of reverse exchange, 
only important in hetero epitaxy. The surface diffusion 
through exchange diminishes the probability of pure Co 
islands formation, which are the responsible of the mul- 
tiple height and vacancy island formation. In the same 
way, the exchange mechanism suppresses the formation 
of the binding alloyed steps, responsible of the step deco- 
ration. Additionally, the coverage of the Co by Pb atoms 
diminishes the effect of the different surface energy be- 
tween Cu and Co stabilizing even the pure Co island. 
This effect turns irrelevant the non zero probability of 
the formation of pure Co island. 
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FIG. 2. Co comparative island stability over Cu(lll) with 
and without a surfactant layer. Upper (first column) and lat- 
eral (second column) view. Upper panel: Original situation 
for a Co island over the critical size. Middle panel: snapshot 
of the Co island over nude Cu, the multilayer island forma- 
tion, together with the pools of vacancies can be observed. 
Lower panel: the evolution of the same Co island, under the 
surfactant layers is shown. The island, as well as the sub- 
strate, shows some kind of corrugation, but for comparative 
times, the island stability is proved 



FIG. 1. Surface diffusion of Cu and Co over Cu(lll) 
pre-covered with a monolayer of Pb obtained by means of MC. 
First column, x, y and z evolution of Cu over Pb/Cu(lll). 
The first drop corresponds to the initial Cu Pb exchange; Sec- 
ond column: as in a, for Co adatom diffusion. Lower panel: 
top views showing the exchange between the Cu substrate 
atoms (gray clear - gray obscuro circle) and the subsurface 
evolving adatom, i.e. Cu or Co (black circle). Note that 
while the exchange Cu/Cu is with the next-nearest neighbor, 
in the Co case is with the nearest one. The Pb atoms are not 
shown in the last frames for clarity 
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